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Abstract 



(1.1) 



Convergence of Rothe's method for the fully nonlinear parabolic 
equation ut + F(D 2 u, Du, u, x, t) = is considered under some continu- 
ity assumptions on F. We show that the Rothe solutions are Lipschitz in 
time, and they solve the equation in the viscosity sense. As an immedi- 
ate corollary we get Lipschitz behavior in time of the viscosity solutions 
CN) . of our equation. 

§ : 1 Introduction 

o 

r-| ■ We consider the boundary value problem: 

u t + F(D l x u, D x u, u, x, t) = in D 

> 

We assume that u(x,t) satisfies Equation ( II .ip in the viscosity sense, that 
D = Q, x [0, T], and that F is uniformly elliptic. We recall that a function u is 
a viscosity subsolution (supers olution) of Equation fl.lj) if for any p 6 C 2,1 (D) 
and any e > which satisfies 

(ft + F(D 2 p, Dip, p, x, t) > e > (12) 
(or (ft + F(D 2 p, Dp, p, x, t) < — e < for the supersolution case) 

u — p cannot attain a local maximum (minimum) of 0. (In other words, p 
cannot touch u from above (below).) u is a solution if it is both a subsolution 
and a supersolution. Uniform ellipticity means that there are positive numbers 
A and A such that for any positive definite matrix iV we have 

-A||iV|| < F(M + N, P, v, x, t)-F(M, P, v, x, t) < -X\\N\ \ , (1.3) 
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or what is equivalent: 



F(M + N, P, v, x, t) < F(M, P, v, x, t) + A||i\T|| - A||iV + || (1.4) 

for all M and N. So the negative Laplacian would fit our structure conditions, 
and not the positive Laplacian. 

For linear F'a Rothe's method is commonly used as a numerical approxima- 
tion. Rothe's method corresponds to doing a backward Euler approximation 
in Banach space, and is also known sometimes as the method of lines. Kikuchi 
and Kacur have studied the convergence properties of Rothe's method in a 
series of recent papers in a variety of function spaces. (See jKj and |KKj and 
the references therein.) Pluschke has studied the quasilinear case (see [P]), 
but the current work appears to be the first extension of Rothe's method to 
viscosity solutions in the fully nonlinear setting. 

We fix our mesh size, < h < 1, and define z n ^{x) recursively. We take 
zo,h(x) = 0, and for n > we let z n+ i t h be the solution of the elliptic problem: 

F(D 2 z n+1A , Dz n+1>h , z n+1>h , x, (n + l)h) + 

%n,h 

~ h 

Zn+i,h(x) = on dQ . t 

z n ,h( x ) is an approximation to u(x, nh), and we define the linear interpolation, 

U h (x,t), by 

\(m + l)h — t] . . \t — mh] . . , 

U h (x, t) := ^ -± l -z m . h (x) + 1 . ^ Z m+l , h {x) , (1.6) 

h h 

where m G IN is chosen so that t G [mh, (m + l)h). When U(x, t) : = 
lim^ Uh(x, t) exists, we call it the Rothe limit. 

We make the following continuity assumption, which corresponds to the 
structure condition (SC) of CCKSj: We assume that for every R > 0, there 
is a nondecreasing continuous function ctr such that uir(0) = 0, and if X and 
Y are symmetric matrices and |r|, \s\ < R, then 

V -(M - N) - 7 |P - Q\ - lo r ((s - r)+) 

< F(M, P, r, x, t) - F(N, Q, s, x, t) (1.7) 

< V + {M — N) + 7 |P -Q\+ uo R {{r - s)+) . 

(We will define the Pucci operators V + and V~ in a moment.) Lastly we 
assume that for any compact subset of H nxn x ]R n x 1R we have 

\\F(M, P, r, ■, •)|| c o^<C<oo. (1.8) 



in Q 



(1.5) 
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2 Elliptic Preliminaries 

To study the elliptic problem in each time we introduce the Pucci extremal 
operators which we give by 

V+(M) = -ATr(M+) + ATr(M-) and 

V (M) = — ATr(M+) + ATr(M-) ^ ' 

where Tr(M) is the trace of M. 

2.1 Lemma (First Bound). There exists a constant C which is independent 
of h such that 

\\zi, h \\ L °°<C . (2.2) 

Proof. By symmetry it suffices to show that there exists a constant C such 
that Z\ t h{x) < C. zx ; h satisfies 

= F(D 2 z hh , Dz 1>h , z 1>h , x, h) + 

h 

> F(D 2 z ljh , Dz lh , z l h , x, h) on {z l h > 0} . 



and therefore by Lemma 2.11 of [CCKSj Z\^ is also a solution of 

V-{D 2 z lth ) - 7|£te ltft | + F(0, 0, z 1)h , x, h) < (2.3) 

on the set where it is positive. By Proposition 2.12 of |CCKSj we now have 

zi, h < C\\(-F(0, 0, z 1>h , x, h)) + \\ L u (r+{zih)) , (2.4) 

where r + denotes the upper contact set. Now by using our structure condi- 
tions, we see that 

-M(-^) + ) < F(0, 0, z lth , x, h) - F(0, 0, 0, x, h) < a R (( Zl , h ) + ) . 

Let f(x) := F(0, 0, 0, x, h). Then we see that 

[-F(0, 0, z 1>h , x, h)) + X{Zhh>0} < \f(x)\ • (2.5) 

By combining Equations ()2.4|) . (j2.5j) . and f)1.8j) we are done. 
Q.E.D. 

2.2 Theorem (First Lipschitz Bound). There exists a constant C which 
is independent of h such that 

||*i,fclU~<C7i. (2.6) 
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Proof. Let the absolute maximum of z\^ be attained at xq, and let 7 : = 

Zi,h{.%o)- We have 

7 = z 1>h (x ) = -hF(D 2 z hh (x ), 0, 7, x , h) 

in the viscosity sense. Since the plane U(x) = 7 touches Z\^ from above at xq 
we know that 

7 < -/iF(0, 0, 7, x , h) (2.7) 

which we know is bounded from above by a constant times h by using the 

previous lemma along with Equations (jl.7|) and (jl.8|) . 

Q.E.D. 



2.3 Theorem (Pucci Inequalities). Let f(x) := F(0, 0, 0, x, (n + l)h). The 

following inequalities are satisfied in the viscosity sense. 



V~(D 2 z n+ i) - jWDzn+xW - a R ((-z n+l )+) + ^±1 + f < 5t 
< V + (D 2 z n+1 ) + j\\Dz n+1 \\ + + Zj f- + / 



(2.8) 



Proof. Let ip touch z n+ \ from above at xq. Since 

F(D 2 z n+1 ,Dz n+ i,z n+1 ,x,(n + l)h) + 

h n 

in the viscosity sense, by using Equation ()1.7|) we have 
j- > F(D 2 ip, Dip, ip, x, (n + l)h) + I 

> p-(DV) - 711^11 " M("V?) + ) + \ + f ■ 

The other inequality is proven in the same fashion by touching z n+ \ from be- 
low. 
Q.E.D. 

Now we introduce the sup-convolution z e - of Zj and the inf- convolution zj >e 
of Zj which are defined by 

(\x — y\ 2 \ f \x — y\ 2 \ 

z o(y) — ) Zj ' e ^ := lei ( Zj ^ + — 2e — J ^'^ 

The next two lemmas are taken from [J] and |JLSj respectively: 

2.4 Theorem (Basic Sup-Convolution Properties). Let C IR n be 

bounded, let u G C{Q), and let u e and u € be its sup and inf- convolution re- 
spectively. Then 
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(i) u e ,u e G C ' 1 ^). 

(ii) M e J, m and u e t u as e [ uniformly on Q . 

(Hi) For every e > 0, there are measurable functions M e , M t : Q — > «S(n) snc/i 

n £ (|/) = M £ (i) + <D« e (i), y-x>+i<M e (x)(y-x), y-x>+o(|a;-?/| 2 ) 

for a.e. x E Q (and similarly with u e and M e ) where S{n) is the set of 
n x n symmetric matrices. 

(iv) M e (x) > -(l/e)I and M t (x) < (l/e)I for a.e x G SI. 

(v) Ifu v ' £ and are standard mollifications of u e and u e respectively, then 
D 2 u r >' e > -(l/e)J and D 2 u^ -> M e {x) a.e. x G ft as r] -> 0, and 
similarly with D 2 u^. 

2.5 Theorem (Sup-Convolutions are Subsolutions). If u is a bounded 
viscosity subsolution of 

F(D 2 u, Du, u, x) = f(x) in B 1 

and f, F are continuous, then 

F(M e (x), Du e (x), u(x e ), x e ) < f(x e ) a.e. in B^^^i/* 

where x e G B\ is any point which satisfies 

I cc^ tic I 

u\x) = u{x e ) . 

Obviously inf-convolutions are supersolutions. 

2.6 Remark. What we have denoted by u x e " is typically denoted by "x*," 
but we need to emphasize the dependence on e to avoid confusion later. (For 
inf-convolutions we will use u x e ") 

2.7 Remark. Note that 

I r^>€ I 2 

= u £ (x e ) - u £ (x) < 2| |u| |oo , 

so that \x t — x\ < 2(e| |m| loo) 1 / 2 . In particular, if x G B 1 _ 2 ( e \\ u \\ oc )i/2 then x £ G B ± , 
explaining the appearance of this set above. 
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By using Equation ()1.5|) and the previous theorem, we get 
F(M^ +lth (x), Dz* n+1>h (x), x\ (n + l)h) + ^ 



Z n+l,h\ X 



< 



(2-10) 



h 



for a.e. x. Similarly we have 



F(M ni h, e (x), Dz nthte (x), Z nAt (x e ), X e , Till) + 



> 



h 

^n— l,h (-^e) 

h 



(2.11) 



Henceforth we suppress the h in the subscripts. By taking a difference we get: 

F(M< +1 (x), Dz* n+1 (x), z^ n+1 (x*),x% (n + l)/i) 



- F(M n>e (x), Dz nie (x), z n7 e(x e ), x e ,nh) + 



< 



h 



(2.12) 



Now we express the difference clS cL telescoping sum. 



F(M^ +1 (x), ^(i), < +1 (x e ), x% (n + l)/i) 

- F(M nie (x), Dz n! e(x), Z ntt (x t ),X t , nil) 



= F (M* +1 (x) , Dzl +1 (x) , < +1 (x* ) , x* 
- F(M n , e (x),D< +1 (x),< +1 (^),x e 

+ F (M„ je (a;) , Dz^ +1 (x) , z^ +1 (x e ) , x e 
- F(M n , e (x),Dz 

+ F(M n>e (x),Dz n>e {x),z e n+1 (x e ),x* 
- F(M nte (x),Dz ntt (x),z nte (x e ),x e 



n + l)h) 
n + l)h) 

n + l)h) 
n + l)h) 

n + l)h) 
n + l)h) 



(2.13) 



+ F(M n!e (x), Dz n>e (x), z nte (x e ),x e , (n + l)h) 
- F(M n;e (x), Dz n , £ (x), z n!e (x e ),x e , nh) 



For the time being we will assume that F is differentiate in the variables we 
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need in order to define 



<HM := Jo F Xij (tM^ +1 (x) + (1 -t)M n , e (x), 

Dz e n+1 (x), < +1 (x e ), x e , (n + l)h) dt 

bi(x) : = Jl F Pt (M n , e (x),tDz* n+1 (x) + (1 - t)Dz n , e {x) 
z f n+1 (x e ), x e , (n + l)h) dt 



(2.14) 



c[x) 



(1 — t)z n>e (x e ), x e , (n + 



Now we use the fundamental theorem on each pair of terms to get the following 
equalities: 



F(M^ +1 (x),Dzl +1 (x),z^ n+1 (x^,x^(n + l)h) ) 
- F{M n>e {x), Dz e n+1 (x), < +1 (x e ), x\ (n + l)h) 



a iA X ) [ M n+Uij)( X ) ~ M n,e,(ij)( X ) 



(2.15) 



F(M n>e (x), Dz* n+1 (x), z* n+1 (x% x\ (n + l)h) 

- F(M nje (x), Dz n>e (x), < +1 (x e ), x\ (n + l)h) 

= bi(x)Di (z e n+x {x) - Z n>e (x)) 

F(M nte (x), Dz n>e (x), z e n+1 (x e ),x e , (n + l)h) 

- F(M nie (x), Dz n ,e(x), Z ntt (x t ),X e , (n + l)fl) 



C(X 



n+1 



X 



Zn,e(.%e)) • 



By the continuity assumption (jl.8j) we also know 

\F(M n>e (x), Dz n>e (x), z Hje (x e ), x e , (n + l)h) 
-F(M n>e (x), Dz n>e (x), z Hje (x e ), x e , nh)\ 



< 



<J2{h) 



(2.16) 



(2.17) 



(2.18) 



where the <7j are moduli of continuity. By defining c(x) := c(x) + (1/h), by 
defining w e n (x) := z* +1 (x) - z n>€ (x) , by defining D jjW e n : = M^ +1 - M n>e (x), and 
by using the equations beginning with Equation (J2.12j) we can conclude 



a ij (x)D ij w e n (x) + bi(x)DiW e n (x) + c(x) [w e n (x e ) + z n , e (x e ) - z ntt (x e )\ 

Z n [X ) Z n -i\X e j , , , . 

< h oAe) + a 2 {h) tor a.e. x 

h 



(2.19) 
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We define a semiconvex function to be a function which becomes convex if 
7||x|| 2 is added to the function, and 7 is sufficiently large. In particular, by 
Theorem (12 .4j) (iv) we know that sup-convolutions are semiconvex, as is the 
function w e n in Equation (I2.19J) . Now we need the following results (see p. 56 
Theorem A. 2 and p. 58 Lemma A. 3 of [CILJ). 



2.8 Theorem (Aleksandrov's Theorem). Semiconvex functions are twice 
differentiable a.e. 

2.9 Lemma (Jensen's Lemma). Let if : H n — > 1R be semiconvex, and let 
x be a strict local maximum off. For p G ]R n , set <f p {x) = <p{x)+ < p,x > . 
Then for any r, 5 > 0, the set 

K := {x G B r {x) : there exists p G B$ such that ip p has a local max at x} 

has positive measure. 

We also state for future use the following estimate: 

2.10 Lemma (Discrete Gronwall Inequality). Assuming {vi}, {B^}, and 
{Di} are sequences of nonnegative numbers which satisfy Vi + % < B(Vi + D+, we 
have 



n—l n—l 

Vn < VP J] B i + ^ 

i=0 i=0 



n-1 



A II B i 



(2.20) 



With these results we can turn to the estimate of the maximum of w n {x) 
z n +i(x) - z n (x). 



3 Rothe Limits Are Lipschitz in Time 

3.1 Theorem (Inductive Estimation). If F satisfies fli. 7| ), then we have 

||w n ||L°°(fi) < lkn-i|U°°(n) + ha 2 (h) . (3.1) 

Furthermore, if we assume that cr 2 (/i) < Ch, then for any fixed T > 0, we have 
the Lipschitz estimate: 

< C , (3.2) 

L°°(n) 

where C is independent of h, and [T/h] is the greatest integer less than or 
equal to T/h. 

Proof. We will first prove Equation ( I3.1J1 . To start we will need a few 
assumptions: First that F is differentiable so that Equation ( 12.19)) is valid, 
second that w n has a strict maximum at x, and finally that our approximating 



W[ T /h] 
h 



functions, are twice differentiable at their corresponding maxima which 
we call x e . We will show how to avoid these assumptions after proving the 
estimate in this case. 

Because x is a maximum, and w e n is twice differentiable at x e , Equation 
(PHI!} becomes: 



C(x) [W e n (x e ) + Z n>e (x e ) - Z nje {x e )\ 



Z n \X ) Z n -\\X t ) 

< r H o-i (e) + <72\h) 

h 



(3.3) 



Now by using Remark ()2.7|) we see that 

\x*-x e \ <2(e\\z e n+1 \\ 0O ) 1 / 2 + 2(e\\z n , e \\ 0O ) 1 / 2 , (3.4) 

and then by sending e to zero and using the Holder regularity of the z n in 
space we can conclude: 

c(x)w n (x) + ^ = c(x)w n (x) < + a 2 (h) . (3.5) 

h h 

By using Equation ( ll.7j) we see that F r > and therefore c > 0. So, since 
w n {x) > 0, we have 

< \\w n -x\\ L ao(p) + /ia 2 (/i) . (3.6) 

Obviously we can argue similarly for a negative minimum, but it remains to 
show how to eliminate our regularity assumptions. First of all we need stability 
of our approximations as we take smooth approximators F m , to our original 
F. The desired stability can be found in section 6 of |CILj . (See Remark 
6.3 in |C1L| in particular.) In terms of applying the estimates we conclude 
for our -F m 's for our original F there is no problem, because no matter how 
badly behaved the a^-, bi, and c are in a given approximation to F, they get 
thrown away at the maximum of w n by considerations which depend only on 
the signs of the terms involved. Stated differently, even if the derivatives of 
the approximations are diverging, they are necessarily diverging "in the right 
direction." (Equation f !3.3|) is obtained from Equation ( I2.19j) by using the 
nonpositivity of the second derivative and the vanishing of the gradient of w n 
at its maximum.) 

Next, if w n has a maximum at x but it is not the unique maximum, then 
we look at the equation satisfied by w n (x) := w n (x) — e\x — x\ 2 , where e > 
is arbitrarily small. (Here we need the continuity assumptions on our function 
F given in (|1.7jh ) Finally, to ensure twice differentiability of the w e n at x e we 
appeal to Jensen's Lemma ()2.9|) and Aleksandrov's Theorem 
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Now we need to prove Equation |3.2j) . We apply Lemma Ij2.10j) to Equation 
(EH1) to get 

||w n ||i/x. < ||wo||l°° +nha 2 {h) . (3.7) 

with n = [T/h] and by using Theorem (J2.2j) we are done. 
Q.E.D. 



3.2 Theorem (Existence of the Rothe Limit). If cr 2 (h) < Ch the Rothe 
limit exists and is locally Lipschitz in time. 

Proof. The last theorem shows that the sequence of approximators are 
uniformly Lipschitz in time. In order to show that the approximators are 
equicontinuous in space, we observe that our uniform Lipschitz estimate im- 
plies F(D 2 z n+ i, Dz n+ i, z n+ i, x, (n + l)h) is uniformly bounded in L°° in the 
viscosity sense. Then we invoke the spatial Holder estimates known for fully 
nonlinear elliptic equations (see Theorem 5.1 Part III of [T]). Finally, we can 
invoke Arzela-Ascoli to guarantee the existence of the limit. 
Q.E.D. 



4 Rothe Limits Are Viscosity Solutions 

4.1 Theorem (Rothe Limits Are Viscosity Solutions). The Rothe limit 
is a viscosity solution. 

Proof. By symmetry, it will suffice to prove that the Rothe limit is a viscosity 
subsolution. Suppose not. Then there exists an e > and a supersolution 
if G C 2,l (D) which touches the limit U from above at a point (xo,to) and 
which satisfies: 

(ft + F(D 2 (f, Dip, <p, x, t) > e > , (4.1) 

pointwise in a neighborhood of (xo,to)- By adding a very small multiple of 
\t — t 1 2 + |x — xo\ 2 to (p, by using the continuity of F, and by allowing a 
slightly smaller (but still positive) e in Equation (j4.1j) . we can assume without 
loss of generality that at least for a small neighborhood of (xq, t ) which we 
will call N the only contact between f and U is at (xo,t ). Because of the 
uniform continuity of ip t there exists a neighborhood N of (xo,to) such that 
for a sufficiently small h > and for all (x, t) G N we have the estimate 



<f(x,t)-(p(x,t-h) 

(pt{x,t) 



h 



< — (4.2) 
-10 V ; 
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for all h < h. At this point let M denote a compact set of the form B r (xo) 
x [t — 7, to + 7] which is contained in N R N. Let 

S:=Af\{B r/2 (x ) x (t -7/2, t + 7/2)}. (4.3) 

Because y? > U on <S, there exists a e > such that ip > U + e on 5. By the 
uniform convergence of % to {7, for any 5 G (0, e/3) we can be sure that by 
shrinking h if necessary we have 

<p>U h + e-6 on S while | <p(x , t ) - U h (x ,t )\ < 5/4 . (4.4) 

for all h < h. Choose 5 sufficiently small to guarantee that 

MS) < e /2 , (4.5) 
where ujr is the modulus given in Equation 1)1.7)) and e is taken from Equation 

Fix h < mm{h, 7/IO}, and sufficiently small to guarantee that 

\<p(x ,t)-U h {x ,t)\<6/2 (4.6) 

for any t E [t — h,t ]. (Here we need the second half of Equation (14.4)1 and 
uniform continuity in time of <p(x ,t) — Uh(x ,t) as h I 0. The continuity is 
guaranteed by the previous theorem.) If we let if = if + c, for any c with 
|c| < 5, then by using Equations (jl. 7)1 and ()4.5)) we will have 

ft + F(£>V, ^, x, t) > e - uo R (6) > e/2 > . (4.7) 

Now we choose c so that the minimum oiif — Uh on the set jVn{(a;, t) :t G INh} 
is equal to zero. In other words, y5 is allowed to be less than between the 
mesh values of t, but when considering only mesh values of t, we can say that 
if touches Uh from above. 




x 
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Because of Equation ( I4.6|) , we can be sure that |c| < 5/2. Let (xo,mh) be 
a point where (p = Uh- By using Equation ( 14 .4|) we know that (xo,mh) G 
5 r/2 (^o) x (t - 7 /2, t + 7/2). 

If we return to our z„ we now have the following setting: 

(a) tp(x,mh) touches z m (x) from above at x~q. 

(b) 0(x, (m — l)h) > z m _i(x) in B r / 2 (x ). (Because h < 7/IO we know that 
B r / 2 (x ) x {t = (m - C A/".) 

Because (x~o,mh) C N we have (see Equation ()4.2|) ) 

</5(xo, — ^(^o, (m — 1)/;,) 



<pt(x ,mh) 



< — (4.8) 

-10 V ; 



h 

Recall that viscosity solution of 

* m + F ( D2 ^> x > mh)=0. (4.9) 

Since (p(x,mh) touches z m (x) from above at x , we conclude that 

y ~^' l(Xo) + F(D 2 (p, Dtp, <p, £0, mh)<0. (4.10) 

Now we combine Equations (|4.7|) , (jUTUj) , and (|4.8jl to get 

e/2 < fi t (x ) + F(D 2 (p,D(p,(p,x ,mh) 
z m -i(x ) - <p(x , mh) 



< <pt(x ) + 

< <ft{x ) + 



<^(fo, (m — l)h) — <p(x~o, mh) 



h 

< e/10 



which is a contradiction. 
Q.E.D. 



4.2 Corollary (Lipschitz Regularity in Time). There is a solution to 
Equation which is Lipschitz in time. 

Proof. Simply combine the two previous results. 
Q.E.D. 
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